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Abstract 

We present a covariant formalism for studying nonlinear perturba- 
tions of scalar fields. In particular, we consider the case of two scalar 
fields and introduce the notion of adiabatic and isocurvature covec- 
tors. We obtain differential equations governing the evolution of these 
two covectors, as well as the evolution equation for the covector asso- 
ciated with the curvature perturbation. The form of these equations 
is very close to the analogous equations obtained in the linear theory, 
but our equations are fully nonlinear and exact. As an application 
of our formalism, we expand these equations at second order in the 
perturbations. On large scales, we obtain a closed system of coupled 
scalar equations giving the evolution of the second-order adiabatic and 
entropy perturbations in terms of the first-order perturbations. These 
equations in general contain a nonlocal term which, however, decays 
rapidly in an expanding universe. 
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1 Introduction 



In a recent series of articles [1, 2, 3] we have developed a new formalism for 
dealing with nonlinear cosmological perturbations in a covariant approach. 
As illustrated in our previous works, the main advantages of this new ap- 
proach are its (relative) simplicity and its broad range of validity. Indeed, 
most of the earlier studies of nonlinear cosmological perturbations rely on 
some approximation from the start, either by considering quantities only up 
to second order in a perturbative approach [4, 5, 6, 7, 8, 9, 10, 11, 12] or by 
restricting their application to super-Hubble scales [13, 14, 15, 16, 17, 18, 19, 
20, 21, 22], or both [23]. By contrast, our approach relies on exact equations. 

Furthermore, by defining perturbations in a covariant and geometrical 
way, our approach avoids many of the technical difficulties encountered in 
the previous approaches, such as the problem of the physical interpreta- 
tion of second-order perturbations in the usual coordinate-based perturba- 
tion theory. Because we do not make any approximation, it is also rather 
straightforward to make the connection with the other approaches. 

In our previous papers we have focused our attention on fluids, either 
perfect or dissipative. The purpose of the present work is to extend our 
formalism to scalar fields, as they play an essential role in early universe 
models. 

An elegant covariant approach to treat perfect fluid inhomogeneities about 
a Friedmann-Lemaitre- Robertson- Walker (FLRW) spacetime was introduced 
in [24], based on earlier work by Hawking [25]. This approach was later 
extended to scalar field dominated universes in [26], although essentially by 
describing the scalar field as an effective fluid. 

Here we consider the situation where several scalar fields are present, 
concentrating on the case of two scalar fields. In the latter case, it has been 
shown in the linear limit that one can define so-called adiabatic and entropy 
components for the perturbations by a suitable rotation of the original scalar 
field perturbations [27, 28]. These two components obey a coupled system 
of second-order (in time) differential equations, which decouple when the 
background evolution in field space corresponds to a straight trajectory. This 
decomposition can also be useful for analyzing the generation of adiabatic and 
entropy perturbations during inflation, which in general can be correlated, 
as first pointed out in [29]. 

A nonlinear formalism for scalar fields with a similar decomposition into 
adiabatic and entropy components has been developed in [20], making use of 
the long wavelength approximation, which consists in dropping higher-order 
terms in the spatial gradients, assuming that these are negligible on large 
scales [13, 14, 15]. In the present work, we show that it is possible to gener- 
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alize this decomposition into adiabatic and isocurvature components to the 
exact and fully nonlinear situation, without resorting to any approximation. 
In order to do so, we introduce adiabatic and isocurvature covectors, which 
generalize the linear definitions. We then derive a coupled system of second- 
order (with respect to time) equations for these two components which are 
exact and covariant, and mimic the equations of the linear theory. 

The covariant equations that we obtain can be expanded at first and sec- 
ond order in the perturbations by working in some coordinate system. At 
first order in the perturbations, it is straightforward to recover the results 
of [27], i.e., two coupled evolution equations for the adiabatic and entropy 
perturbations and an evolution equation for the curvature perturbation on 
comoving or uniform density hypersurfaces. At second order in the pertur- 
bations, we extend the results of [11, 30] and we derive, on large scales, an 
evolution equation for the (second-order) entropy perturbation, which sources 
the evolution of the (second-order) curvature perturbation. 

Furthermore, we will show that the second-order two-field system has 
qualitatively new features with respect to the linear case or the second-order 
single-field case. In particular, the scalar equations that we derive contain 
generically a nonlocal term due to the impossibility of writing the total mo- 
mentum of the two fields as a total gradient. However, we will show that 
this nonlocal term can be neglected, on large scales, due to the expansion of 
the universe. 

This work is organized as follows: in the next section we present the basic 
equations to treat a system of several scalar fields in the covariant formalism. 
In Sec. 3 we discuss the perturbations of a single-field dominated universe, 
while in Sec. 4 we extend this treatment to two fields. In Sec. 5 we give a 
set of covariant equations that describe the two-field system in two limits: 
the linear regime and the large-scale regime. From our covariant equations 
we derive, in the coordinate approach, the evolution equations for the linear 
perturbations in Sec. 6 and for the second-order perturbations on large scales 
in Sec. 7. Finally, in Sec. 8 we draw our conclusions. 

2 Covariant formalism for several scalar fields 

Let us consider an arbitrary unit timelike vector u a = dx a /dr {u a u a = — 1), 
which defines a congruence of cosmological observers. The spatial projection 
tensor orthogonal to the four-velocity u a is defined by 

hab = g a b + u a u b , {h a b h b c = h a c , h a b u b = 0). (1) 
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To describe the time evolution, our covariant definition of the time derivative 
throughout this work will be the Lie derivative with respect to u a , which is 
defined for a generic covector X a by (see e.g. [31]) 

X a EE C u X a EE u b V b X a + X b V a u\ (2) 

and will be denoted by a dot. We emphasize that this notation differs from 
the convention used in the literature based on the covariant formalism (see 
e.g. [24]), as well as in our previous papers [1, 2, 3], where a dot denoted 
simply the covariant derivation along u a . The two definitions coincide only 
for scalar quantities, 

/ = u b V b f. (3) 

To describe perturbations in the covariant approach, we introduce the 
projected covariant derivative orthogonal to the four-velocity u a , D a , in the 
spirit of [24]. For a generic tensor, the definition is 

D a T b c : =h a % e ...h f c ...W d T e f r . (4) 

In particular, for a scalar quantity /, the definition reduces to 

D a f = hy b f = d a f + uJ. (5) 

We also introduce the familiar decomposition 

VfeW a = o- ab + u ab + ^Qh ab - a a u b , (6) 

with the (symmetric) shear tensor a ab and the (antisymmetric) vorticity ten- 
sor uj ab \ the volume expansion, 9, is defined by 

6 EE V a U a , (7) 

while 

a a = u b V b u a (8) 

is the acceleration vector. 

Let us now consider iV scalar fields minimally coupled to gravity with 
Lagrangian density 

C = -\d a <p I dP<p I -V(<p K ), (9) 

where V is the potential. We assume, for simplicity, canonical kinetic terms. 
Here and in the following, summation over the field indices (I, J, . . .) will be 
implicit. The energy- momentum tensor derived from this Lagrangian reads 

T ab = datpfiw 1 - l -g ab (Uy^y + 2V) . (10) 
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Given an arbitrary unit timelike vector field u a , it is always possible to de- 
compose the total energy momentum tensor as 

Tab = {p + P)u a u b + q a u b + u a q b + g ab P + ir ab , (11) 

where p, P, q a and n ab are respectively the energy density, pressure, mo- 
mentum and anisotropic stress tensor measured in the frame defined by u a . 
Starting from the energy-momentum tensor (10) one finds 

p = T ab u a u b = l - (W + D a ^D a Vl ) + V, (12) 

P = l -h ac T ab h b c = l - (W - ^A^V) - V, (13) 
q a = -u b T bc h c a = -tprDaip 1 , (14) 
vr ab = KT cd h d b - Ph ab = D aVl D b ^ - hi ab D cVl D c ^. (15) 

The evolution equations for the scalar fields are obtained from the varia- 
tion of the action with respect to the scalar fields. One thus gets N Klein- 
Gordon equations, given by 

dV 

V„VVj = j—. (16) 

d(fi 

Like for the energy-momentum tensor, it is useful to consider a decomposition 
into (covariant) time-like and space-like gradients defined with respect to u a . 
One then finds that the above Klein-Gordon equations can be re-expressed 
in the form 

^ j + Qifj + V m - D a D a Vl - a a D a ifj = 0, (17) 
where we have defined 

dV 

V n , w . (18) 

In the following sections, we will consider in turn the case of a single scalar 
field, which has already been studied briefly in our covariant approach [1], 
and the case of two scalar fields. 



3 Single field 

3.1 Equation of motion 

For a single scalar field, which will be denoted by 0, the Klein-Gordon equa- 
tion (17) reads 

dV 

(f> + ®<P+-T--D a D a ( i ) -a a D a( j) = 0. (19) 
a<p 



Although this equation is fully nonlinear, it closely mimics the homogeneous 
Klein-Gordon equation in a FLRW spacetime when one identifies 0/3 with 
the local Hubble parameter. However, the last two terms depend explicitly on 
the spatial gradients but it is possible, in the case of a single-field, to get rid 
of them by a particular choice of the unit vector u a . Indeed, if the gradient 
V a is time-like, which is the case if the spacetime is sufficiently close to 
FLRW, then one can always choose for u a the unit vector field orthogonal to 
the hypersurfaces, 

< om = ± V (20) 

i-V b <pV h (t) 



where the sign of the right hand side depends whether V a is future or past 
oriented. In this comoving frame, the spatial gradient of 0, i.e. D a (p, auto- 
matically vanishes. Therefore, for this particular choice of u a , the last two 
terms on the left hand side of (19) disappear and the Klein-Gordon equation 
looks formally similar to the homogeneous version although it remains fully 
inhomogeneous and nonlinear. 

One can go one step further by deriving an exact and covariant equation 
that mimics the equation of motion governing the linear perturbations of a 
scalar field in a perturbed FLRW spacetime. The trick is to consider the 
gradient of the scalar field, i.e., the covector 

0a = V a 0. (21) 

This can be decomposed into a spatial gradient and a longitudinal compo- 
nent, 

a = D a <j) - u a . (22) 

Note that with the choice (20) the spatial gradient disappears in the above ex- 
pression and this equation can be used to rewrite the comoving four-velocity 
as < om = -V a 0/0. 

We now construct a second-order (in time) evolution equation for <f> a . We 
recall that a dot stands for the Lie derivative along u a , as defined in Eq. (2). 
We can derive the evolution equation for a by taking the spacetime gradient 
of Eq. (19) and noting that, for any scalar 0, the Lie derivative with respect 
to u a and the spacetime gradient (but not the spatial gradient) commute, 
i.e., 

V a 0=(V a 0)\ (23) 

Equation (19) then yields 

0« + O0a + V H <p a = -0V a O + V a {D b D b <t)) + V a (a 6 A,0). (24) 

This equation mimics the analogous perturbation equation at linear order but 
also incorporates the fully nonlinear dynamics of the scalar field perturbation. 
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3.2 Integrated expansion perturbation on comoving sli- 
ces 

In [1, 2] it was shown that one can define a covariant generalization of the 
comoving curvature perturbation by using appropriate combinations of spa- 
tially projected gradients. For a scalar field, a natural choice is the covariant 
integrated expansion perturbation on comoving hypersurfaces lZ a , defined as 

a 

U a = -D a a + -jD a <f>, (25) 
where a is the integrated volume expansion along u a , 

a = - fdrQ (0 = 3d). (26) 
3 J 

Since 0/3 corresponds to the local Hubble parameter, one sees that the 
quantity a can be interpreted as the number of e-folds measured along the 
world-line of a cosmological observer with four-velocity u a . 

When the four-velocity u a is chosen to be comoving with the scalar field, 
as defined in Eq. (20), then the last term in the definition (25) drops out. For 
this particular case, the evolution equation for 1Z a has already been given in 
Ref. [1, 2]. 

Note a useful property of TZ a : one can replace in its definition (25) the 
spatial gradients D a by partial or covariant derivatives, 

a 

U a = -V a a + T a . (27) 

<P 

The same property applies to the nonlinear generalization of the comoving 
Sasaki- Mukhanov variable for a scalar field, which can be defined as 

Q a = D a ct>-^D a a = ^n a . (28) 
a a 

3.3 Integrated expansion perturbation on uniform en- 
ergy density slices 

It is also possible to generalize the curvature perturbation on uniform energy 
density hypersurfaces, as shown in detail in our previous works [1, 2]. The 
key role is then played by the covector £ a defined as 

( a = D a a - -D aP . (29) 
P 



7 



For a perfect fluid, the quantity ( a satisfies a remarkably simple first-order 
evolution equation [f, 2], 

e 2 

(a = jrT a , (30) 

(the dot stands, as before, for a Lie derivative with respect to u a ) where 

T a = D a P-jD aP (31) 

is the nonlinear nonadiabatic pressure perturbation. For a barotropic fluid, 
T a = and ( a is conserved on all scales. The relation (30) for ( a can be seen 
as a generalization of the familiar conservation law for (, the linear curvature 
perturbation on uniform energy hypersurfaces. It can also be extended to a 
non-perfect fluid [3]. 

Note however that the relation between our ( a and the traditional ( is 
somewhat subtle on small scales. As we showed explicitly in [2], if we work in 
some coordinate system the spatial components of ( a correspond essentially 
to the spatial gradient of ( on super-Hubble scales, but on smaller scales there 
is a non-negligible extra term. Consequently, on small scales Eq. (30) does 
not imply, for adiabatic perturbations, the conservation of the traditional £. 

For a scalar field, the comoving and uniform density integrated expansion 
perturbations ( a and lZ a satisfy 

(a + n a = ~ I^D aP - tD a <^ , (32) 

which simply follows from their respective definitions. The right hand side 
can be interpreted as the "shift" between hypersurfaces of constant p and 
hypersurfaces of constant <fi and the term inside the parenthesis represents 
the nonlinear generalization of the so-called comoving energy density of a 
single scalar field. 

By choosing u a = u® om as defined in Eq. (20), the energy-momentum 
tensor of a single scalar field can be written in the perfect fluid form, i.e., 
with vanishing q a and n a b and the energy density p and pressure P given by 

P = \tf + V, P = \tf - V, (33) 

as can be checked by specializing Eqs. (12-15) to a single field and setting 
»«0 = 0. 

The scalar field cannot, strictly speaking, be considered as a barotropic 
fluid. Indeed, Eq. (33) implies 

P = p-2V, (34) 
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which, after substitution in the definition (31) of the nonadiabatic pressure 
co vector r a , gives 




where we have used D a V = 0. Substituting the time derivative of Eq. (34), 
one gets 

r a = 2|^IV, (36) 
and the evolution equation of ( a for a single scalar field is then given by 

C. = l^D a p, (37) 

where we have used p = — @0 2 to obtain this expression. As we will show in 
Sec. 5.2, the right hand side of this equation can be neglected on large scales, 
so that ( a is conserved in this limit. 



4 Two scalar fields 

We now consider several scalar fields and for simplicity we restrict our anal- 
ysis to the case of two scalar fields, which we will denote by and \- 



4.1 Adiabatic and entropy covectors 

In the two-field case, it is possible to introduce a particular basis in the 
field space in which various field dependent quantities are decomposed into 
so-called adiabatic and entropy components. In the linear theory, this de- 
composition was first introduced in [27] for two fields. For the multi-field 
case, it is discussed in [28] in the linear theory and in [20] in the nonlinear 
context. 

In our case, the corresponding basis consists, in the two-dimensional field 
space, of a unit vector e£ defined in the direction of the velocity of the two 
fields, and thus tangent to the trajectory in field space, and of a unit vector 
el defined along the direction orthogonal to it, namely 

e[ = 1 (0, *) , el = 1 (-*, $) . (38) 
yj 2 + X 2 \I4> 2 + X 2 

An immediate consequence of the above definitions is the identity 

8 j = e I a e aJ + ele sJ , (39) 



9 



which will be very useful in the derivation of many relations. 

To make the notation shorter it is convenient to introduce the angle 9 
defined by 

cos# = ^=i_, sm9 = —?^=, (40) 
V <P 2 + X 2 V 4> 2 + X 2 



so that 



= (cos#,sin6>) , = (- sin 9, cos 9) . (41) 

This angle, in contrast with the linear theory case where it is a background 
quantity that depends only on time, is here an inhomogeneous quantity which 
depends on time and space. By taking the time derivative of the basis vectors 
and ef, we get 

ei = 0el, el = -H- ( 42 ) 
It is also convenient to introduce the formal notation 



& = 4¥+I 2 - (43) 

Note however that this notation can be misleading as, in general, in the 
nonlinear context, a is not the derivative along u a of a scalar field o. 

Making use of the basis (38), one can then introduce two linear combina- 
tions of the scalar field gradients and thus define two covectors, respectively 
denoted by a a and s a , as 

a a ee ei Va^fi = cos 9 V a + sin 9 V a X, (44) 
s a = e^V a v?/ = -sin#V a + cos#V a x. (45) 

We will call these two covectors the adiabatic and entropy covectors, re- 
spectively, by analogy with the similar definitions in the linear context [27]. 
Whereas the entropy covector s a is orthogonal to the four- velocity u a , i.e., 
u a s a = 0, this is not the case for a a which contains a "longitudinal" compo- 
nent: u a a a = a. It is also useful to introduce the spatially projected version 
of (44-45), 

a a = e i D aVl = °a + 0~U a , EE e^D^I = S a . (46) 

Let us now consider the "adiabatic" combination of the Klein-Gordon 
equations, i.e., the contraction of (17) by e. l a . By noting that 

e&i = cr, e-i^i = (47) 

and by defining 

VV = eX, (48) 
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one obtains 

a + B& + V a - eiD a D a Vl - a a a^ = 0. (49) 
The fourth term can be rewritten as 

eiD a D a Vl = D a (e^Vj) - (D a ei)D a Vl = D a a^ - {e sI D a e I a )s\ (50) 

where we have used the definition of a^, Eq. (46), and the identity (39). By 
noting that 

e sI V a e T a = -e aI V a el = -(s a + 9a a ), (51) 

o 

one finds 

elD a D a Vl = D a a--Y (s) , (52) 

where we have defined 

Y {s) = ^(s a + 9ai)s a . (53) 

Thus, finally, the adiabatic combination of the Klein-Gordon equations can 
be written as 

ij + e& + V a = V a <J^-Y (s) , (54) 
where we have used the property 

DX+flV^VX, (55) 

which is valid for any covector orthogonal to u a . 

Let us now consider the "entropic" combination of the Klein-Gordon 
equations, i.e., the contraction of (17) with e\. By using 

e{(pi = Oct, (56) 

and by defining the entropic gradient of the potential, 

V s = e T s V m , (57) 

one finds 

&9 + V s - eiDaD^j - a a s a = 0. (58) 
One can rewrite the last two terms by using the identity 

elD a D a Vl = D a s a + Y {a) , (59) 

with 

Y {a) = ^(s a + 9ai)a^ : (60) 
11 



and by applying the property (55) to the covector s a . One finally gets 

&9 + V s = V a s a + Y (a) . (61) 

To summarize, we have been able to replace the Klein-Gordon equations 
for the fields <fi and \ by two new equations which embody the evolution 
along the adiabatic and entropy directions respectively. The left hand side 
of these equations has exactly the same form as the homogeneous equations 
in a Friedmann universe. Although our covariant equations (54) and (61) 
look very similar to the homogeneous equations, they capture the fully non- 
linear dynamics of the scalar fields. This is manifest in the right hand side 
of these equations, which contain nonlinear (quadratic) terms represented 
by F( s ) and Y"( CT ), sourcing the adiabatic and entropy equations respectively. 
Whereas these equations generalize the background evolution equations, we 
will go one step further in the next subsection by deriving fully nonlinear 
and exact equations which mimic and generalize the linearized equations for 
the adiabatic and entropy components. 



4.2 Evolution of the adiabatic and entropy covectors 

We now derive evolution equations for the covectors o~ a and s a . More pre- 
cisely, our purpose is to find two evolution equations, which are second order 
in time (with respect to the Lie derivative along u a ) and which mimic the 
equations obtained in the linear theory for the perturbations 5a and 8s (see 
[27]). 

Let us begin with the evolution equation for o~ a . Starting from its def- 
inition (44), one finds that its time derivative, i.e., the Lie derivative with 
respect to u a , is given by 

& a = eiV a ¥i + 0s a = V a a + 9s a , (62) 

where the last equality is obtained by using (pj = ae aI . A further time 
derivative yields 

o a = V a a + 9s a + 9s a . (63) 

The next step consists in using (54) to eliminate a in the above expression. 
This gives 

a a + && a + dv a e + (v >aa + e^) °a - v a (vv) = (e~Y) 

+ (6 - V as + 60) s a - V a Y (s) , (64) 
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where we have used the relation 

V a V a = V aa a a + V as s a + ^(s a + 9a a ), (65) 

a 

and introduced the notation 

V = e T e J V V = e I e J V V = e I e J V (66) 

for the second derivatives of the potential. 

The evolution equation for a a can be decomposed into a longitudinal part, 
obtained by contracting (64) with u a , and an orthogonal part obtained by 
contraction with h a b- By using the relation 

V* = V >7IJ a + eV s , (67) 

it is not difficult to see that the longitudinal part yields in fact the time 
derivative of (54). What is more interesting is the orthogonal or spatial part 
which can be written as 

(a^ + 0(cr a ) ± + aD a e + (V :(Ta + °l ~ D a (V c ^) = *a 

+ (e - v as + ee) s a - D a Y {s) . (68) 

Let us now consider the evolution equation for s a . From Eq. (51) the 
time derivative of s a is given by 

s a = &elV a e aI - Qo a . (69) 

Taking another time derivative, one finds 

s a = ~0a a - 9a a + -{s a + 6a a ) + &V a 0, (70) 
a 

where we have used (51) and (42). We can now use the entropic equation 
(61) to get rid of V a 9. Furthermore, using the relation 

V a V s = V sa a a + V ss s a - ^(s a + 9a a ), (71) 

a 

we finally obtain 

-s a --(a + V a )s a + (V ss - 9 2 )s a - V a (V c s c ) = -29& a 
a 



+ 



-{a + V a )-9-V 
a 



o a + V a F (CT) . (72) 
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As for the adiabatic equation, the longitudinal part of this equation, upon 
using the relation 

V s = V as a - 9V a , (73) 
yields the time derivative of Eq. (61). The orthogonal part, instead, yields 

s a - \{<J + V a )s a + (V ss - 9 2 )s a - D a (V c s c ) = -26{a a ) L 

+ D a Y {a) , (74) 

where we have used the property that the covectors s a and s a are purely 
spatial, i.e., that (sa)" 1 = s a and (sa)" 1 = s a . 

Starting from the fully nonlinear Klein-Gordon equations, we have thus 
managed to obtain a system of two coupled equations (68) and (74), which 
govern the evolution of our nonlinear adiabatic and entropy components. 
These equations are one among the main results of the present work. Re- 
markably, they are rather simple and they look very similar to their linearized 
counterparts. As we will see later, immediately deduced from these equations 
are the already known evolution equations for the linear adiabatic and en- 
tropy components. Furthermore, since our equations are exact, they can 
be used to go beyond the linear order, up to second or higher orders. To 
illustrate this, we will show explicitly in Sec. 7 how to derive from these 
equations, in a systematic way, the evolution of the second-order adiabatic 
and entropy components. 




4.3 Generalized covariant perturbations 

In this subsection, we will be first interested in the covariant generalization 
of the comoving energy density and curvature perturbations in the context of 
a two-field system. We will then consider the generalization of the curvature 
perturbation on uniform energy density hypersurfaces. 
Let us introduce the covector 

e a = D a p-^, (75) 
a 

which can be interpreted as a covariant generalization of the comoving energy 
density perturbation. In order to obtain the explicit expression of e a in terms 
of a a and s a , let us rewrite the components (12-15) of the energy-momentum 
tensor in the form 

p=i(a 2 + n)+\/ (76) 
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P = I _ l n ) - V, (77) 

q a = -&<t£, (78) 

7r ab = Il ab - -h ab U, (79) 

where we have defined 

Tl ab = a^at + s a s b , U = a^a^ c + s c s c . (80) 
Substituting the expression for p into the definition (75), one gets 

e a = ~ + (V s - 6&)s a + l - (^D a U - ^a^j , (81) 

where we have used the relation 

D a & = cr a + <yu a - 9s a , (82) 

which is a consequence of (62). 

As noticed above, the longitudinal projection of Eq. (72) yields the time 
derivative of (61), which reads 

-V a -6- V as = -a- -(D c s c + Y {a) )\ (83) 

This relation, together with Eqs. (54) and (81), can be employed to rewrite 
Eq. (74) in the form 



s a + 



" \ (V c ^ - Y( S ))] s a + (v ySS + 6 2 - 20^) s a - D a (V c s c ) = 
-2^e a + Q - {dJ\ - ^a^j - ^(D c s c + Y {a) ) + D a Y {a) . (84) 

As in the linear theory (see [27]), this gives us an alternative expression for the 
evolution equation for s a , in which the comoving energy density perturbation 
appears explicitly on the right hand side. This expression will be useful in 
Sec. 5.2 when discussing the large-scale evolution of s a . 

Together with the comoving energy density, the comoving curvature per- 
turbation can generalized. For the general case of several scalar fields, this 
is done by defining the comoving integrated expansion perturbation, 

a 

Tia = ~D a a - J^rj^q a - (85) 
15 



The definition of the Sasaki- Mukhanov covector given for a single scalar field 
in Eq. (28) can then be extended to the case of several fields, by defining for 
each field 

Qi = D& 1 - ^D a a. (86) 
a 

Thus, the comoving covector lZ a can also be written as 

K " = T^¥) (87) 

In the two- field case, the definition (85) reduces, using (78), to 

n a = -D a a + -a^. (88) 
a 

Furthermore, one can generalize the Sasaki-Mukhanov variable to the adia- 
batic covector by defining 

Q a = e aI Q T a = a£- -D a a. (89) 
a 

Instead of giving the evolution equation of lZ a we will resort to a fluid 
description by considering the covariant generalization of the uniform den- 
sity curvature perturbation, i.e., the integrated expansion perturbation on 
uniform density hypersurfaces ( a , defined in Eq. (29). In the two-field case, 
these two quantities are related by 

Ca + n a = —e a . (90) 
P 

In contrast with the case of a single scalar field, which can always be 
described as a perfect fluid, the total energy-momentum for two or more 
scalar fields will in general correspond to that of a dissipative fluid and the 
nonlinear formalism developed in [3] will thus be useful in this case. 

The adiabatic Klein-Gordon equation (54) can be rewritten as a conti- 
nuity equation for the total energy density (76) and pressure (77), which 
reads 

p + q(p + p) = v, (91) 

with the "dissipative" term 

v = & (vv^ - y {s) ) + ien + In. (92) 

In [3] it was shown that the evolution equation for ( a for a dissipative 
fluid, which generalizes (30), is given by 

(a = |? (r a + S a ) , (93) 

3p 
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where the second source term on the right hand side, due to the dissipative 
nature of the fluid, is defined in terms of V as 



S a = -i (r> a V - jD^ + ^ (# a e - ®D a p^j . (94) 

We have called this term the dissipative nonadiabatic pressure perturbation 
because it has the same form as Y a , provided one replaces the pressure P by 
a dissipative pressure —V/Q defined in terms of dissipative fluid quantities. 

We now concentrate on the covector T a and rewrite it by taking into 
account the fact that our effective fluid consists of two scalar fields. In this 
case, the pressure and energy density are related by 

P = p-2V- 2 -U ) (95) 

as can be inferred from Eqs. (76-77). Inserting this relation into the definition 
of T a , Eq. (31), one gets 

r a = ^1 - D aP - 2D a V - ^D a U. (96) 

Using the time derivative of (95) and the relation 

D a V = V„at + V s s a , (97) 

one finally obtains 

T a = 2jV a e a - 2V s s a - 1 (^D a U - j D a p^j , (98) 

where we have introduced e a defined in Eq. (75). The above equation ex- 
presses the nonlinear nonadiabatic pressure perturbation in the two-field case, 
which sources Eq. (93). It will be useful in the next section, where we will 
consider both the linear and super-Hubble approximations of our evolution 
equations. 

The above equation expresses the nonlinear nonadiabatic pressure pertur- 
bation in the two-field case. The evolution of £„, governed by Eq. (93), is thus 
sourced by a rather complicated term obtained by summing Y a in Eq. (98) 
and E a in Eq. (94). As we will see in the next section, this term simplifies 
considerably when either the linear or the super-Hubble approximations are 
taken. 
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5 Approximate equations 



In this section we study the evolution equations of the two-field system under 
two types of approximations: the linear limit and the limit where we neglect 
higher orders in spatial gradients. In an expanding FLRW universe, the 
latter corresponds to the large scale limit. In the rest of the paper we will 
use the symbol ~ to denote an equality at the linear level, and the symbol 
~ to denote an equality valid only on large scales. 

5.1 Homogeneous and linearized equations 

In many cosmological applications, since our Universe appears to be close 
to a FLRW universe on large scales, it is sufficient to restrict oneself to the 
linearized version of the evolution equations. Before working within some 
generic coordinate system in the next section, we first consider this lineariza- 
tion procedure directly at the level of the covariant equations, as in [24]. 
In a strictly FLRW universe, all the spatial gradients vanish and therefore 

a£ = 0, s a = 0. (FLRW) (99) 

Consequently, the scalar quantities F( s ) and Yjv), defined respectively in (53) 
and (60), vanish and the evolution equations for a and 9, respectively (54) 
and (61), reduce to 

a + 3H& + V a = 0, (FLRW) (100) 

and 

ab + V s = 0, (FLRW) (101) 

where we have introduced the Hubble parameter H = 0/3. Not surprisingly, 
the above equations correspond exactly to the homogeneous equations given 
in [27]. Furthermore, in a FLRW universe, all the terms in Eq. (74) for s a 
vanish. This is the same for Eq. (68) for a a projected orthogonally to u a . 
Although its longitudinal component does not vanish it is simply the time 
derivative of (100). 

At linearized order, we treat the covectors and s a , which vanish at 
zeroth order, as first-order quantities. Similarly, their derivatives s a , s a , 
(6-a) 1 - and ((To) 1 - are first-order quantities. Therefore the linearized version of 
the evolution equations are simply obtained by keeping only the homogeneous 
terms in the coefficients multiplying the spatial projection of o~ a , s a and their 
derivatives. The linearized versions of (68) and (84) are thus, respectively, 

(aa) 1 - + 3H(& a ) L + &D a e + (V aa - 9 2 ) ai - D a (D c a^) ~ 2 (9s a ) ' - 2^6s a , 

V02) 
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and 

s a + 3Hs a + (V ss + 39 2 ) s a - D a (D c s c ) ~ -2-e a . (103) 

Note that the terms involving Y( s ), Yjy) and II have disappeared, since these 
scalars are quadratic in first-order quantities. We have also replaced V c o^ 
by D c cr^, as well as V c s c by D c s c , since their difference is quadratic in first- 
order quantities, according to (55). Indeed, the acceleration vector a b , which 
vanishes at zeroth order, is considered as a first-order quantity. Finally, we 
have used Eqs. (100-101) (together with the time derivative of the latter 
equation), to rewrite the right hand side of the linearized equations in the 
above form. 

One can also linearize the evolution equation for ( a . As discussed above, 
the terms containing Y^ and II can be neglected. The dissipative term V 
thus reduces to 

V ~ &D a a^, (104) 
while the expression (98) for T a becomes simply 

T a ~ 2jV a e a - 2V s s a . (105) 

Therefore, the evolution equation for ( a , Eq. (93), can be written, at linear 
order, as 

L ^ ^V a e a + ^V s s a + ±rD a {D b ai). (106) 

As expected, the linearized equations (102-103) and (106) are equivalent 
to the linear equations obtained in the coordinate based approach in [27]. 
This will be even more explicit in Sec. 6 where we introduce a generic coor- 
dinate system and compute explicitly the components of our various tensors 
and equations. 

5.2 Expansion in spatial gradients 

Apart from the linearization procedure, there is another approximation in the 
cosmological context which applies to describe the Universe on very large 
scales, even beyond the present Hubble radius where in principle it could 
strongly deviate from a FLRW universe. This approximation is based on an 
expansion in spatial gradients, which are small for scales larger than the local 
Hubble radius [13, 14, 15]. 

In this perspective one sees, from their definition (46), that and s a are 
first-order quantities with respect to spatial gradients because they are linear 
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combinations of spatial gradients. The scalars F( s ) and are second order 
with respect to spatial gradients since they are quadratic in and s a (or 
their time derivatives). Hence, the right hand side of Eq. (54) and of Eq. (61) 
can be neglected on large scales, so that these two equations become, in the 
large-scale limit, 

a + && + V a « 0, (107) 

and 

9^-^. (108) 
a 

Although they look very similar to the homogeneous equations (100) and 
(101), these equations are fully inhomogeneous and encode the evolution of 
nonlinearities on large scales. This limit illustrates the separate universe 
picture [32, 33] where the inhomogeneous universe can be described, on large 
scales, as juxtaposed Friedmann homogeneous universes. 

If, so far, the order in spatial gradients seems to coincide with the per- 
turbative classification of the previous subsection, it differs however for the 
term 

V^=Z}V + ^, (109) 

which is first order perturbatively but second order in spatial gradients, at 
least for the first term on the right hand side, since is already first order 
in spatial gradients. For the second term, we will assume that u a can be 
chosen so that a c is at least first order in spatial gradients. We will show this 
explicitly in the next section by working in a coordinate system. 

With these prescriptions, the evolution equation of and s a obtained 
at lowest order in spatial gradients become 

(a a ) ± + e^) 1 + &D a e + (v aa - e 2 ) ^ » 2(e Sa ) - 2e^f Sa , (no) 

and 

S a + Q-Sa + (V ss + 36 2 ) S a W -2~6 a , (111) 

where we have dropped the terms containing II in (84), which are two orders 
higher than s a in spatial gradients. 

We now expand the evolution equation for ( a , Eq. (93), by neglecting 
higher-order spatial gradients in the two terms on the right hand side of this 
equation. The nonadiabatic pressure perturbation becomes 

r a w 2-V a e a - 2V s s a , (112) 
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while the dissipative nonadiabatic pressure perturbation S a can be com- 
pletely dropped, since the dissipative term V is at least second order in the 
spatial gradients and thus S a is third order in the spatial gradients. Equation 
(93) therefore becomes, on large scales, 



2 26 
(a~ T 7 1 V a e a + —V s s a . (113) 



Note that the lowest order limit in spatial gradients of the evolution equations 
of (Ja, s a and ( a , respectively Eqs. (110-111) and (113), are similar to their 
linear counterparts except for the terms D a (D c a^) and D a (D c s c ), which are 
third order in spatial gradients and therefore negligible in the spatial gradient 
expansion. This is because in these equations the terms that are higher than 
linear order in the perturbative expansion turn out to be also higher than 
first order in the spatial gradient expansion. 

We now concentrate our attention on the comoving energy density per- 
turbation, e a , defined in Eq. (75). Until now we have only made use of the 
Klein-Gordon equations of the scalar fields. However, in order to study the 
behavior of the comoving energy density, we will now make use of the Einstein 
equations, in particular of the so-called constraint equations. The projection 
of the Einstein equations along u a yields the energy constraint in a covariant 
form, 

u a G ab u b = 8irGp. (114) 

If one assumes that the vector field u a is hypersurface orthogonal, it is possi- 
ble to use the Gauss-Codacci equations (see e.g., [31]) and the decomposition 

D b u a = a ab + - Oh ab , (115) 

which is the spatially projected version of (6) (with uj ab = since u a is here 
hypersurface orthogonal), in order to rewrite the energy constraint as 

1 f ^R + -6 2 - o- ab a ab ) = 8tiG P , (116) 



2 V 3 

where is the intrinsic Ricci scalar of the space-like hypersurfaces orthog- 
onal to u a . 

The mixed projection of Einstein's equations yields the covariant momen- 
tum constraint 

u b G bc h c a = 8nGq a , (117) 
which can be rewritten, via Gauss-Codacci relations and Eq. (115), as 

D b a a b -^D a Q = 8nGq a . (118) 
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By combining the energy and momentum constraints, one obtains the non- 
linear covariant version of the generalized Poisson equation, which in the 
linear theory relates the comoving energy density to the Bardeen's potential 
defined from the curvature perturbation. Here one finds 

l -D a - a bc a bc ) + QD b a b a = 8nGe a , (119) 

where we have introduced on the right hand side 

e a = D a p - Bq a = D a p + Gaa±. (120) 

This quantity can be seen as an alternative generalization of the comoving 
energy density perturbation, since, in the linear limit, it is equivalent to e a 
defined in (75). In the fully nonlinear case, the two quantities differ and, 
using (91), one finds 

e a -e a = ^(v- ien) (121) 

This difference becomes however negligible on large scales, 

e a « e . (122) 

Now, the left hand side of Eq. (119) contains the projected gradient of the 
Ricci scalar, D a ^R. From its definition in terms of derivatives of the metric 
[31], it can be shown that, for a perturbed FLRW universe, this term is of 
third order in the spatial gradients. Equation (119) can thus be used to show 
that e a in Eqs. (90), (111) and (113) can be neglected on large scales, if the 
shear can also be neglected in this limit. Indeed, on large scales, the shear 
rapidly decreases in an expanding perturbed FLRW universe. Thus, in this 
limit the comoving and uniform density integrated expansion perturbations 
( a and 1Z a coincide (up to a sign), 

Ca + n a w 0. (123) 

Furthermore, one can rewrite Eqs. (Ill) and (113) as a closed coupled system 
of equations, describing the large-scale nonlinear evolution of adiabatic and 
entropy perturbations, 

's a + ®s a + (V ss + W 2 ) s a w 0, (124) 

and 

2 G 2 

Ca » — V>a- (125) 

3 p 

These two equations are analogous to the equations derived for nonlinear 
perturbations by Rigopoulos et al. [21] in the long wavelength approximation. 
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6 Linear perturbations 



We now relate our covariant approach with the more familiar coordinate 
based formalism. We first examine the linear perturbations in the present 
section and we consider second-order perturbations in the next one. 

Let us thus introduce generic coordinates x M = {t,x 1 } to describe an 
almost-FLRW spacetime. Here a prime will denote a partial derivative with 
respect to the cosmic time t, i.e. ' = d/dt, since the dot has been reserved 
to denote the Lie derivative with respect to u a . 

The background spacetime is a FLRW spacetime, endowed with the met- 
ric 

ds 2 = g^dx^dx" = -dt 2 + a(t) 2 % j dx i dx j . (126) 
At linear order, the spacetime geometry is described by the perturbed metric 

ds 2 = {g„ v + Sg^) dx^dx\ (127) 

where the components of the metric perturbations can be written as 

^oo = -2A, 5g 0i = aBi, Sgij = a 2 H i:j . (128) 

As usual we decompose Bi and Hij in the forms 

B t = V t B + BY, (129) 
= -2^ ltJ + 2ViVjE + 2V {l El + 2£j, (130) 

where Bj and are transverse, i.e., WiB Vl = = WiE v \ and Efj is 
transverse and traceless, i.e., WiE T%: > = and — 0- Here V* denotes 

the three-dimensional covariant derivative with respect to the homogeneous 
spatial metric 7^' (which is also used to lower or raise the spatial indices). The 
matter fields are similarly decomposed into a background and a perturbed 
part, 

<p I (t,x i )=<p I (t) + 6<p I (t,x i ). (131) 

We now need to specify the components of the unit vector u a , which 
defines the time derivation in our covariant approach. At zeroth order, it is 
of course natural to take it orthogonal to the homogeneous slices. At first 
order we choose, for simplicity, w M such that Ui — 0. This implies that, up to 
first order, the components of are given by 

u" = {1 - A, -B^a}, (132) 

and those of the "acceleration" vector are given by 

^ = {0,VVL/a 2 } . (133) 
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This confirms that d 1 can be considered as first order in spatial gradients, in 
agreement with our assumption of the previous section. 

Since our formalism relies on many covectors, it is useful to first consider 
a generic covector X a and work out the components of its time derivative X a . 
To make the explicit calculation, it is convenient to replace in the definition 
of the Lie derivative (2) the covariant derivatives by partial derivatives and 
write 

X a = u b d b X a + X b d a u b . (134) 
At zeroth order, the components of X a are simply 

X^{X' , 0,0,0}, (135) 

assuming that the spatial components Xi vanish so as to respect the symme- 
tries of the geometry. At first order, we get from (134) and (132) 

5(X ) = 5X' - (X A)', 8(Xi) = 5X[ - X d t A. (136) 

These results can be applied to the particular case of the covector a = 
V a 0, which we introduced earlier. Its components are, by definition, 

O = 0', 0, = 0, (137) 

at zeroth order and 

50o = 50', Hi = dM, (138) 
at first order. Specializing the relations (136) to X a = a , one gets 

5(0 O ) = 50" - 4>'A' - fi'A, 6(fa) = d t (50' - 4>'A) . (139) 

Applying once more (136), now with X a = a , one finds 

6(4>i) = di (50" - 4>'A' - 20"A) . (140) 

We will not need the time component of a , which involves third order deriva- 
tives with respect to cosmic time. 

By using the above results one finds that, at linear order, the spatial 
components of Eq. (24) for a single scalar field correspond to the spatial 
gradient of the following scalar equation: 



50" + 3#50' + V M 8(f) - ^V 2 50 = -2V^A + 0' \A' + 3^' - V 2 (£' - B/a) 

a 



;i4i) 
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where we have used 

e = 3H, 5& = -3HA-3ip' + V 2 (E'-B/a), (142) 

as well as the background equation 

4>" + 3H4>' + V t4> = . (143) 

In the case of several scalar fields, the same equation applies for each </?/ with 
the replacement of V^Sc/) by V tLpiV>J 5ip J . 

Let us now consider, for the special case of two scalar fields, the adiabatic 
and entropic covectors a a and s a , which we have introduced earlier. Note that 
with the choice of four-velocity (132) aj~ = Oi. The background equations of 
motion can be deduced immediately from Eqs. (54) and (61) and read 

a" + 3Ha' + V a = 0, (144) 

a'B' + V s = 0, (145) 

with 

a = a' = v/0' 2 + x' 2 . (146) 

From its definition, Eq. (44), one finds that the spatial components of a a at 
linear order can be expressed as 

6ai = ^.dM + %di8 X = dM, (147) 
a a 

with the notation 

5a=^5<f>+^5 X , (148) 
a a 

which coincides with the definition of [27]. By using Eq. (136) with (146-147) 
one finds that the spatial components of the first and second time derivatives 
are given by 

6(&i) = di {5a' - a' A) , 6(<7i) = d t (5a" - a' A' - 2a" A) . (149) 
The same procedure for s a gives 

5si = di5s, 5s = — 5x — —.54>, (150) 
<j a 

which also coincides with the notation of [27]. Since s a , in contrast to a a , 
has no longitudinal component, So = and the spatial components of s a and 
s a are simply 

6(s i ) = d i 6s', 6(s i ) = d i 6s". (151) 
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Plugging the explicit components (149) and (151) into the linearized equa- 
tions for a a and s a , given by (102) and (103) respectively, one obtains easily 
the linearized equations for 5a and 5s. These read, respectively, 



5a" + 3H5a' 



6' 2 ) 5a - \f 2 5a = 2(6'5s)' - 2^&5s 



and 



-2V a A + a' A' + 3^' - V 2 (E' - B/a) 



5s" + 3H5s' + (V ss + 36' 2 )5s - -^V 2 5s = -2—Se . 

a 2 a' 



(152) 



(153) 



In the latter equation we have introduced the first-order comoving energy 
density perturbation 5e, defined by 



5ei = di5e, 5e = 5p — —5a, 

a' 

which follows from the definition (75) of e a . Using 

K>2 



(154) 



(155) 



5p = a' 5a - Aa' 2 + VJa + 2V s 5s, 

one sees that 5e can be expressed as 

5e = a 5a' - a 2 A - a" 5a + 2V s 5s . (156) 

Moreover, linearizing the spatial components of the energy constraint 
(116) yields 



3H (if/ + HA) - —V 2 ip + H(a E' - aB) = -4ttG Sp, 



(157) 



while the momentum constraint (118) gives, since 5qt = —a'diSa, 

iP' + HA = AnGa'Sa. (158) 

Combining the two above constraints yields the relativistic Poisson-like equa- 
tion 



1 



ip + H(a E — aB) = AttG 5e, 



(159) 



which can also be directly obtained by linearizing the spatial components of 
Eq. (119). This equation shows that the comoving energy density perturba- 
tion 5e is second order in the spatial gradients, and thus negligible on large 
scales in Eq. (153). 



26 



In contrast with 5s, the quantity 5a is not gauge-invariant. This is why 
it is useful to consider the gauge invariant Sasaki-Mukhanov variable Qsm, 
defined as [34] 

Q SM = 5a + (160) 



Note that the above traditional definition does not follow exactly from our 
definition of Q a given earlier in Eq. (89). Indeed, from Q a , one can extract 
a scalar quantity Q defined as 



Qi diQ , 



a 



Q = 5a - —5a, 
ti 



(161) 



where 5a can be written in terms of metric perturbations by making use of 
Eqs. (26) and (142) (see [2]), 



Sa = -if; + - I V 2 (E' - B/a)dt. 



(162) 



Thus the scalar variable Q coincides with Qsm only in the large-scale limit. 
In the flat gauge, defined by ip = 0, 5a coincides with Qsm, 



5a = Q 



SM- 



(163) 



In this gauge, it is possible to use the momentum constraint equations (157) 
to derive the metric perturbation A as a function of Qsm, 



A = 



H' 



Ha 



-,QsM, 



(164) 



and one can write the Poisson equation (159) as 



V 2 (E'-B/a) 



H' 



Ha' 



Qsm ~ ^0'5s 



(165) 



where we have used the expression (156) specialized to the flat gauge. By 
replacing Eqs. (164) and (159) into the evolution equation of 5a, one finds 
the evolution equation of Qsm [35, 27], 



Qsm + SHQ'sm + 



e' 2 -r—\^- + ~—--\- — 



H' fV a H 



a" 



H \ a 



H a' 



= 2(8' 6s)' -29'[^- + 



IT 
H 



5s. 



Qsm 
(166) 
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When one considers only large scales, the expression (165) reduces to 

Q'sm <2sm - 26' '5s » 0, (167) 

which means that there exists a first integral for the quantity Qsm and that 
the second-order equation of motion (166) is not necessary in this limit. In 
fact, one can easily check that the large-scale limit of (166) is an automatic 
consequence of the first integral (167). 

Let us now consider the evolution equation for ( a with two scalar fields, 
Eq. (106). The spatial components of ( a , at linear order, are given by [2] 

5( t = diC, C = 5a- ^5p. (168) 

P 

From Eq. (162) the scalar variable ( is thus related to the Bardeen gauge 
invariant variable (b, defined as [36, 37] 

Cb = - %&P, (169) 
P 

in a similar way to how Q is related to Qsm, and on large scales these two 
quantities coincide. 

According to Eq. (136), the spatial components of ( a are 

S(Ci) = d^'. (170) 

One then finds that the spatial components of Eq. (106) correspond to the 
spatial gradient of 

C = ^ffc - ™?8s + ^r~ 2 ^5a. (171) 
6a' A a' 6<y'a l 

On large scales, the spatial gradients can be neglected and the above equation 
reduces to the well-known relation 

C w - 2 jr^' 5s - ( 172 ) 

Eq. (171) can be easily specialized to the case of a single scalar field by simply 
setting 8a = 5(f) and 5s = 0. One recovers in particular that on large scales ( 
is conserved. Note also that from Eq. (90) there is a simple relation between 
C and 71, 

( + K = -^5e, (173) 

which shows that ( and —1Z coincide on large scales. 

In conclusion, all the equations in this section, derived directly from our 
covariant formalism, exactly reproduce the linear results of [27]. In the next 
section we turn to the second-order perturbations and derive novel results. 
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7 Second order perturbations 



In this section we consider second-order perturbations. We will thus decom- 
pose any scalar quantity X as 

X(t,x*) = X(t) + SXWfax*) + SXWfax*), (174) 

where X(t) is the background part and SX^ and 5X^ are respectively the 
first and second-order contributions (note that we do not follow here the 
convention of including a numerical factor 1/2 in front of the second-order 
contribution). In our subsequent equations, to simplify the notation, we will 
often omit the index W for the first-order quantities, unless it is required for 
clarity reasons. 

Our main purpose will be to expand our equations governing a a , s a and 
( a at second order in the perturbations. Before undertaking this task, it is 
instructive to recall the results of our previous works, in particular of [2] , on 
the second-order component of ( a for a general fluid. For a discussion on the 
expansion of Q and the gauge invariance of ( at order higher than second see 
[38]. 



7.1 ( a at second order and the issue of gauge- invariance 

As shown in [2], the second-order expression of the spatial components Q can 
be written, after some manipulations, in the form 



cf = *c (2) + ^c (1)/ 

p 



with 



P P 



+ 2 



y) s " 



(175) 



(176) 



and being given in Eq. (168). In the large-scale limit, i.e., keeping only 
the scalar perturbations without gradients, where 



a ~ In a — tp — ip 2 , 



(177) 



we have also shown that our second-order quantity could be easily related 
to the second-order quantity defined by Malik and Wands [9], which we 
denote here by Cmw> namer y 



a(2) c^s _ /- 
s ~ ww s: 



(2) 



.(1) ^ 
MW 



(178) 
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Since £mw was constructed by explicitly requiring gauge-invariance, this im- 
plies that, on large scales, our also behaves like a gauge-invariant quan- 
tity. It is instructive to understand directly why our is indeed gauge- 
invariant on large scales. 

Under a second-order coordinate transformation, 

- 5" = x» - % + \^, v - $)> (179) 

generated by the vector fields £fo and £? 2 )j the ^ rs ^ an d second-order pertur- 
bations of a tensor T transform as [4] 

ST (i) 5T«+^ (1) T(°), - 5T( 2 )+^ (2) T(°) + i4 (i) T(°)+£ e(1) 5T«. 

(180) 

Since £ a vanishes at zeroth order, ( a is automatically gauge-invariant at first 
order, according to the first expression above. However, ( a is not gauge- 
invariant at second order and the corresponding gauge transformation is 
given, according to Eq. (180), by 

+ (181) 

Concentrating now on large scales, one finds from the expression of the Lie 
derivative that 

A (1) C^^i)5oC (1) , (182) 

where we have neglected the terms of higher order in spatial gradients. Con- 
sequently, the second-order spatial components of ( a transform on large scales 
according to 

C (2) - Cf + £(i)<9*C (1) ' (large scales) . (183) 
By noting that, at first order, 

<«*> 

it is easy to see that the quantity 

Cf - ^C (1) ', (185) 
P 

or equivalently c^ 2 -*, is gauge-invariant at second order, on large scales. 
Thus, this proves that defined above is indeed gauge-invariant on large 
scales. 
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7.2 Adiabatic and entropy fields 

Here we derive the evolution equations for the adiabatic and entropy field 
perturbations at second order. For simplicity, we will restrict ourselves to 
large scales and we will thus start from the equations expanded in spatial 
gradients discussed in Sec. 5.2. For convenience, we have collected in the 
appendix various background and first-order expressions that will be used in 
the rest of this section. 

The second-order evolution for the perturbations of a single scalar field, in 
a coordinate based approach, has been considered in [5, 10, 7, 11, 39, 40]. The 
multi-field case, in the large-scale limit, has been studied in detail by Malik 
in [30] (see also [41, 42, 43]) and, using the separate universe approach, in 
[23, 44]. However, the second-order decomposition into adiabatic and entropy 
components has not been given. This decomposition, which appears quite 
involved in the coordinate based approach, becomes natural in our nonlinear 
formalism for scalar fields and can be derived straightforwardly since we have 
already identified the fully nonlinear adiabatic and entropy components. 

We start by expanding the definition of a a in Eq. (44) at second order. 
After some straightforward manipulations, one can write the spatial compo- 
nents of <r,- as 



x ( 2 ) 1 



di(<f)'5<f) {2) + x'$X {2) ) + (Ss' + 9' 8a) d l Ss] . (186) 



To deal with the term Ss'diSs, which cannot be written as a total gradient, 
it is convenient to introduce the spatial vector 

Vi = ^(SsdiSs - Ss'diSs) = -Ss'diSs + ^di(SsSs'), (187) 

which vanishes when Ss' and Ss have the same spatial dependence, i.e., Ss' = 
f(t)Ss. 

By expanding also the definition of s a in Eq. (45), one finds, for Sj and 
<7j, respectively, 

= diSo- {2) + —SadiSs - —Vi, (188) 



6v® 

a' 

Ssf = diSs^ + 5 4-diSs', (189) 
a 

with 

5a m = ^0(2) + ^ Sx (2) + A-ssSs', (190) 
Ss m = _^ (2) + l Sx (2) _ ^ ( 5s > + 9 ' 5a \ . (191) 
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The form of the right hand side of Eq. (189) has been chosen by analogy with 
the form (175). Since s a vanishes at zeroth order, arguments similar to those 
of the previous subsection ensure that 5s {2 \ defined in Eq. (191), is gauge 
invariant on large scales. The form of 5s^ is in some sense dictated by our 
covariant definition. 

Note that 5sf^ contains the first-order adiabatic perturbation. This is due 
to the fact that the adiabatic and entropy components are defined locally: 
whereas the first-order components are defined with respect to a background 
basis in field space, which is only time dependent, the second-order compo- 
nents will be sensitive to the first-order fluctuations of the field space basis, 
which can be expressed in terms of the first-order adiabatic and entropy com- 
ponents. The adiabatic component a a does not vanish at zeroth order and 
is not a gauge invariant variable. Our formalism does not dictate its 
form, which is thus chosen at our convenience. 

Let us now discuss further the adiabatic component. Since 5a^ is not 
gauge-invariant on large scales, in contrast with 5s^ 2 \ it is useful to consider 
our generalization of the Sasaki-Mukhanov variable, Q a , defined in Eq. (89). 
Its spatial components can be expanded at second order in the perturbations, 
similarly to what we have done with ( a , 



Qf > = d t Q^ + ^Q' + -08,5s - l-Vi, (192) 
where is defined as 



H H 



Q {1) ' + l{^j'sa-0'5s 



(193) 



From this expression and Eq. (177) it is natural to define 



Q^-\[^]^-0'5s 



(194) 



as the local part of the scalar gauge invariant second-order Sasaki-Mukhanov 
variable. Restricted to a single scalar field, this definition coincides with the 
one given in [30]. Note that we cannot write Eq. (192) in the same form 
as Eq. (175) because the last two terms on the right hand side cannot be 
written as a total spatial gradient. 

The second order (in time) evolution of a a is given by Eq. (110). However, 
on large scales we do not need to compute a second-order differential equation 
because the adiabatic evolution is governed by a first integral, as in the linear 
case. This first integral is obtained directly from the constraint equations and 
it is not necessary to expand (110) at second order in the perturbations. 
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In order to compute this first integral we need the second-order energy 
and momentum constraints, which can be derived by expanding Eqs. (116) 
and (118) and by using 

50 {2) « -HA 2 + 3Aij/ - Gipip'. (195) 

On large scales, one can write the second-order energy constraint equation 

as 



3H 



HA (2) + ^(2), + -Itl- 2A(HA + tf/) 

2 H 



-AnG5p {2) , (196) 



where ^p*- 2 -* is given by 

5p {2) » aW (2)/ - a' 2 ^ 2 ) + V^ (2) + 2K s 5s (2) 

- 2A( ( t , 5 ( j' - a' 2 A - a'9'5s) + A p , (197) 

and A p is a quadratic function of 5a^\ 5s^\ and their first derivatives, given 
by 



A p = -^/((fo' - 20'<5s) - 5s' (-^<Ss + 9' 5a 



+ i(y CTCT - fl' 2 )^ 2 + (V ss + 2# /2 )5s 2 + V s Js5a. (198) 



The second-order momentum constraint equation reads 



ft 



/L4< 2 > + V (2)/ + 2^' - -HA 2 - A(HA + 

(2) 



-4irG8q?\ (199) 



where the second-order momentum <5g- is given, from Eq. (78), by 

5 q {2) = -di ( a'5a (2) + -—5a 2 + 9'5a5s] - -SedM + V t . (200) 
y 2 a J a' 

As already noticed in [20] for large scales, 5q^ cannot be written as a total 
gradient when several scalar fields are present. After neglecting 5e on large 
scales in the above equation, this is manifest because of the presence of V*. 
This implies that, in principle, if Vi does not vanish, one cannot define at 
second-order a comoving gauge, i.e., such that 5q^ = and 5qf^ = 0, in 
contrast with the linear theory or the single-field case. 

However, it is instructive to derive the evolution equation for Vi on large 
scales by using the linear evolution equation for 5s, Eq. (153), neglecting the 
gradient term and 5e at first order. One finds 

V! + SHVi = 0, (201) 
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which implies that, in an expanding universe, Vi will decay like a~ 3 and 
rapidly become negligible even if it is nonzero initially. Consequently, in an 
expanding universe one can in practice ignore Vi on large scales and thus 
define, in an approximate sense, a comoving gauge at second order, which 
coincides with 5a^ = = Sa^. In this approximate comoving gauge, the 
momentum (200) can be written as a total gradient. In the rest of the paper, 
in order to remain as general as possible, we will keep the term Vi. 1 

Similarly to the first-order case, it is possible to combine the energy and 
momentum constraint equations and derive the relativistic Poisson-like equa- 
tion analogous to Eq. (159), which corresponds to the expansion, at second 
order and on large scales, of Eq. (119). By expanding Eq. (120) and using 
(122), one has 

8ef ] « d t 5p {2) - ZHSqf ] - 5Q5q t « 0, (202) 

where the last approximate equality is a consequence of Eqs. (142) and (196) 
and confirms our conclusion of Sec. 5.2 in a covariant context, namely that 
we can neglect e a on large scales. The second-order spatial components of e a 
defined in Eq. (75), can be decomposed as 

5ef ] = dM 2) + ^dM l) ' ~ 3HV, (203) 

<7 



with Se^ defined by 




-i 



a' 



(204) 



It is only when Vi is negligible that the quantity ^e^ 2 -* can be interpreted as 
the comoving energy density at second order. Otherwise, as discussed before, 
the comoving gauge cannot be defined. 

Using the decomposition (203) and the fact that is negligible on large 
scales, Eq. (202) can be written as 

d 2 5e {2) w ZH&Vi. (205) 

When Vi is negligible, and only then, one finds that, like at first order, the 
second-order comoving energy density is negligible on large scales, 

5e {2) w 0. (206) 



1 Note that a similar term (defined for several scalar fields) has appeared in a recent 
paper [45] where it was neglected. Here we have provided the explanation for why it can 
be neglected. 
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Having discussed the general properties of the second-order constraint 
equations, we can now derive the evolution equation of the gauge-invariant 
adiabatic component Qsm- Similarly to what we have shown in the previous 
section for the first-order variables, the simplest way to derive an equation 
satisfied by Qsm is to work in the flat gauge ip^ = = ip( 2 \ where 5a^ 
reduces to <3sm> 



5a 



(2) 



VSM- 



In this gauge, Eqs. (196) and (199) reduce, respectively, to 

(2) 



3H 2 (A& - 2i 2 ) 
Hdi (i (2) - ^A 



-AwGSp 



(2) 



(207) 

(208) 
(209) 



Using the first-order constraint equations Eqs. (164) and (167), these can be 
rewritten as 



A& « -JL 

2V 



Q 



(2), , V a 



SM 



29' 8s™ 



+ 



H' 



Ha 

A 



SM 



Ha 



-,Qsm5s + — f 



(210) 



and 

i(2) 



H' 



Ha' 



(2) 
VSM 



H r 



+ -QsmSs + _ ( _ - 3— ) Q 2 M - -<T 



H 



(211) 

This last equation contains a nonlocal term because we have written the 
momentum constraint as a scalar equation while keeping the second-order 
vector Vi defined in (187). By combining these two relations to get rid of 
one obtains the following first integral for Q^: 



(2)/ 



(H' a"\ . 
+ \H QsM 

1 f ~H'V„ 



(2) - 26'Ss^ 



-2- 



- 3 



a' 



H + 



H a' 
H'\ 



1 

2 



a" 



H 



QsmSs + 



-3 
1 



.IT 

H / 
3H + 



3H 
~H , 



IT 

H 



O 2 



d~ 2 d l Vi. 



(212) 



This equation is the second-order equivalent of (167). We have put together 
on the left hand side all the terms depending on purely second-order quanti- 
ties, and on the right hand side all the terms which are quadratic in first-order 



35 



quantities. As in the first-order case, the entropy perturbation sources the 
evolution of the adiabatic perturbation. The nonlocal term containing Vi 
comes from the momentum constraint and is a new feature with respect to 
the first-order case (or the second-order case for a single scalar field). How- 
ever, as we have discussed earlier, it becomes quickly negligible on large scales 
in an expanding universe, in which case the first integral (212) becomes a 
scalar local equation. 

It is also possible to construct a second-order (in time) differential equa- 

(2) 

tion for Qsm that always looks purely local, even in the presence of a non- 
vanishing Vi. Using the evolution equation of Vi, Eq. (201), this can be done 
by taking an appropriate linear combination of (211) and its time derivative 
so that all nonlocal terms cancel each other. One then obtains a scalar equa- 
tion for a linear combination of A^' and which can be transformed into 

(2) 

a second-order differential equation for Q s ^ by substituting the expression 
for A^ given by the energy constraint. 2 But, even then, the initial con- 
ditions for this second-order differential equation must be compatible with 
Einstein's equations and thus satisfy the nonlocal constraint (212). As a con- 
sistency check, we have also verified that this second-order (in time) equation 
is equivalent to the equation obtained by expanding to second-order (in the 
perturbations) the equation (110) for a a , when combined with the constraints 



We can now derive the evolution equation for the entropy perturbation 
Ss^ 2 \ which will be second order in time. Since we are restricting ourselves 
to large scales, we simply expand the spatial components of Eq. (Ill) up to 
second order. This gives 



where we have neglected the gradient of the comoving energy density at first 
order, 5e^ which according to Eq. (159) is subdominant on large scales. 

To proceed, we need the spatial components of the first and second time 
derivatives of the covectors s a . By using Eq. (134) for Sj at second order, 
and ignoring the higher-order terms in the gradient expansion, one obtains 



2 Notc that a similar second-order equation was obtained in [30] by using the second- 
order constraint from the diagonal space-space components of Einstein's equations. It 
can be shown however that this constraint is redundant with the energy and momentum 
constraints which we are using here and the two procedures are thus essentially equivalent. 



(210-211). 



5(s t ) {2) + 3H5(s z ) (2) + (V ss + 36' 2 ) 5s {2) + 50cW 
+ \SV SS + 6$' 6(6)] di5s « -2-5ef\ 



(213) 



(214) 
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and, by applying once more (134), 

6(3^ « 5sf ] " - A'diSs' - 2Ad l 5s". (215) 

We can then substitute the expressions (189) and (214-215) for the second- 
order entropy component and its derivatives. Using various first-order ex- 
pressions summarized in the appendix, we finally find that the equation (213) 
can be written as the spatial gradient of the following scalar equation 

5s^" + 3H5s^' +(V SS + 36' 2 ) 8s™ » --5s' 2 
-| (»" + ~« V f - SsSs' 

-fo,,,- 5 -V,,-9?pj6 S >-2 e -6 C V\ (216) 

Note that this equation is closed, in the sense that only the entropy field 
perturbation appears: even when 5e^ is not negligible, it can be written in 
terms of 5s and 5s by using Eqs. (205) and (187). Thus, on large scales the 
entropy field evolves independently of the adiabatic components, as in the 
linear theory. 



7.3 Generalized uniform density and comoving curva- 
ture perturbations 

We now derive the large-scale evolution equation for ^ 2 \ expanding at second 
order Eq. (113). In the fluid description, it was shown in [9, 2] that on large 
scales the evolution equation for can be written as 

C (2V ~ ^W 2 ) ^rxtf, (217) 

s p + p p + p 1S1 ' v ; 

where can be read from the second-order decomposition of the quantity 
T a for a fluid, defined in Eq. (31), i.e., 

5TP = + ^r«'. (218) 
P 

In the two-scalar field case considered here we must compare this expan- 
sion with the expression for T a given in the large-scale limit by Eq. (112). 
Expanding this equation at second order, using Eqs. (202) and (189), one 
obtains 

5Tf ) » 2^VJef } - 2V S Uss^ + ^d t 5s') - 25V s d t 5s, (219) 
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where we have used that 5e^ is negligible on large scales. Replacing SV S in 
the last term of Eq. (219) by the expression given in the appendix, one can 
rewrite this equation in the form (218), with 

r( 2 ) « 2yVM 2) - 2V s 5s (2) - V ss 5s 2 + ^-SsSs'. (220) 

One can now rewrite the second-order evolution equation for on large 
scales, Eq. (217), as 



c (2)/ - -4 



20'a'Ss^ - (V ss + AO 12 ) 5s 2 + ^5s5s' - ^Se^ 
v ' ' a' 3Ha' 



(221) 



The last term on the right hand side can be re-expressed in terms of the now 
familiar nonlocal term involving V$, which decays quickly in an expanding 
universe. 

It is also useful to express our results in terms of lZ a . The spatial com- 
ponents of lZ a can be decomposed as 

Tlf = d{RV + %d,{R^' - ILv u (222) 
a' o' A 



with 

H An- r 1 / H\ ' _ H 1 

(223) 



a' a' 



The last term in Eq. (222) comes from the fact that, like e a and in contrast to 
( a , lZ a is defined in terms of the spatial momentum which cannot be expressed 
in general as a pure gradient. When this term can be neglected, and only 
then, VS 2 ^ coincides with the second-order comoving curvature perturbation 
defined in [6, 11]. 

It is easy to derive a first-order (in time) evolution equation for TZ^ 
by noting that and are related on large scales. Indeed, expand- 
ing Eq. (123) at second order using (175) and (222), and neglecting terms 
proportional to 5e^\ one gets 

C (2)+ft(2)« JLfcW. (224) 

3a' 2 

When 6e<® is negligible on large scales, like in an expanding universe where 
we can neglect V^, and VS 2 ^ coincide on large scales as in the single- field 
case [11]. However, this is not true in general in the multi-field case if Vi 
cannot be neglected. 
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From this relation and the evolution equation of (( 2 \ Eq. (221), one can 
find a large-scale evolution equation for , 



H 



29 



■'o>SsW - (V ss + AO' 2 ) 5s 2 + ^5s5s' + ( 1 + 



H' N 
3iP 



<fe< 2 > 



(225) 

The second-order uniform adiabatic field perturbation T^ 2 -* can be related on 
large scales to <5sm> by combining Eqs. (194), (223), and (177). One obtains 



H 

IP 



Qsm - —, {Q'sm - o'<> s ) Qsm - ^ ) <5i M 



(226) 



which can be used, together with the linear first integral (167), to show that 
Eq. (225) is equivalent to the first integral (212). 



7.4 Discussion 

In the two-field system at second order, we have found several qualitatively 
new features that are absent in the linear case or in the second-order single- 
field case, and that can be viewed as purely nonlinear effects. Before con- 
cluding, it is worth discussing these features along with the main results of 
this section. 

Here we have derived coupled equations governing the evolution of the 
adiabatic and entropy perturbations at second order, valid in the large-scale 
limit. The entropy perturbation is described by the second-order field pertur- 
bation 5s^ , which is gauge invariant. The adiabatic perturbation, as usual, 
can be described by one out of three gauge invariant variables. Two choices 
are the Sasaki- Mukhanov variable (?sm> defined in (194), and TZ^ 2 \ defined 
in Eq. (223), the latter being related to the former via the nonlinear equa- 
tion (226). Another choice is the curvature perturbation on uniform energy 
density hypersurfaces, represented, on large scales, by (( 2 \ 

The scalar evolution equations for these three variables, respectively Eqs. 
(212), (225) and (221), carry the same information and share the same char- 
acteristics. They are first order (in time) and they are sourced by the first 
and second-order entropy field. They also contain a nonlocal term which ap- 
pears due to the impossibility of writing the two-field momentum as a total 
gradient. However, we have shown that this nonlocal term decays rapidly in 
an expanding universe. When it can be completely neglected, and only then, 
7^( 2 ) corresponds, on large scales, to the comoving curvature perturbation. 
In this case VS 2 ^ and coincide on large scales. 

The time variation of VS 2 ^ and depends only on the first and second- 
order entropy fields. Thus, the full large-scale evolution of second-order per- 
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turbations in the two-field system is solved by considering the evolution equa- 
tion of either or as well as the evolution equations for the first order 
and second order entropy perturbations, respectively Eqs. (153) and (216). 

8 Conclusions 

In the present work, we have developed a covariant formalism that deals 
with fully nonlinear perturbations in a universe dominated by scalar fields. 
In order to do so, we have introduced a unit vector field u a , which defines 
our time direction. In contrast with the case of a perfect fluid where it is 
natural to define u a as the fluid four-velocity, we have here left u a arbitrary. 
In the case of a single scalar field, it might be convenient to take u a as the 
unit vector field orthogonal to the constant scalar field hypersurfaces, so that 
the total momentum vanishes, q a = 0. However, in the case of several scalar 
fields there is in general no such choice. 

For an arbitrary number of scalar fields, we have shown that it is possible 
to rewrite the fully nonlinear Klein-Gordon equation in a form which mimics 
its homogeneous version. By introducing the gradient of the scalar field, it 
is easy to obtain an equation which is exact and fully nonlinear but at the 
same time closely mimics the linearized Klein-Gordon equation. 

We have paid special attention to the case of two scalar fields, by intro- 
ducing two combinations of the scalar field gradients which we call adiabatic 
and isocurvature covectors because they generalize the definitions introduced 
in the linear theory context. Remarkably, it has been possible to derive for 
these covectors two exact and fully nonlinear second-order differential equa- 
tions, where the time derivative is defined covariantly as the Lie derivative 
with respect to u a . We have also derived the evolution equation of the covari- 
ant variable ( a which generalizes, in the spirit of our formalism, the curvature 
perturbation on uniform energy density hypersurfaces, showing that on large 
scales it is sourced only by the entropy covector. 

We have used these nonlinear equations as a starting point to show that 
results previously obtained in the literature with other approaches can be de- 
rived here in a simple and straightforward way. We have also been able to go 
beyond previous works by computing explicitly the evolution of the second- 
order adiabatic and entropy components on large scales. In particular, using 
the second order energy and momentum constraints from Einstein equations, 
we have derived a first integral of motion satisfied by the second-order adi- 
abatic component, which has the particularity to contain a nonlocal term 
that depends on the first-order entropy perturbation and its time derivative. 
However, this nonlocal term goes rapidly to zero in an expanding universe 
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and the second order adiabatic component is, in this limit, governed by a 
local first-order (in time) evolution equation, sourced by terms depending 
on the second order entropy perturbation as well as, quadratically, on the 
first-order entropy perturbation. Both first and second order entropy per- 
turbations satisfy a second order (in time) evolution equation and the full 
system of equations, valid on large scales, is thus closed. 

A Useful identities in a two-field system 

A.l Background identities 

H' = -AixGa' 2 . (227) 



V' a = <j'V aa + 6%, (228) 

V' s = a'V sc - 9%, (229) 

V; as = a'V aas + 9'(V ss -V aa ), (230) 

V'=a'V ssa -29'V as . (231) 



9' = -^f , (232) 



cr 

0" = -V sa + 9' (Y^-^-\ (233) 
9"> = -a% as + 9" 9> (sv a(7 - V ss - 29' 2 - 2—^- + 3H') . 



(234) 



A. 2 First order identities 



5V a = VvJa + V as 5s - 9' (5s' + 9' 5a), (235) 

5V S = V sa 5a + V ss 5s - ^(5s' + 9'5a), (236) 

5V SS = V ss Ja + V sss 5s - 2^- (5s' + 9' 5a) , (237) 

5V as = V aas 5a + V ass 5s + f % - (5s' + 9'5a). (238) 

V a a 



41 



6(&) = 5a' - 9' 5s - a' A, (239) 
5(0) = -1 [V S A + 5V S + 9' (5a' - 9' 5s)] , (240) 

6(9) = -6V as + 9'SQ +(^- (T -\ 5(9) + 20'^f - 2^^.(241) 

\ a a J a a a 

Acknowledgments: We acknowledge Asko Jokinen, David Lyth, Anu- 
pam Mazumdar, Sami Nurmi, Gerasimos Rigopoulos, Bartjan van Tent and 
David Wands for interesting discussions. FV wishes to thank Ruth Durrer 
and the department of Theoretical Physics of the University of Geneva, David 
Lyth and the department of Physics at Lancaster during the workshop Non- 
Gaussianity from Inflation, 5-9 June 2006, and the Institut d'Astrophysique 
de Paris for their kind hospitality while completing this work. Finally, we 
would like to thank Sebastien Renaux-Petel and Gianmassimo Tasinato for 
pointing out a typo in Eq. (213) and a term missing in Eqs. (220), (221) and 
(225). 



References 

[1] D. Langlois and F. Vernizzi, Phys. Rev. Lett. 95, 091303 (2005) 
[arXiv:astro-ph/0503416] . 

[2] D. Langlois and F. Vernizzi, Phys. Rev. D 72, 103501 (2005) 
[arXiv:astro-ph/0509078] . 

[3] D. Langlois and F. Vernizzi, JCAP 0602, 014 (2006) [arXiv:astro- 
ph/0601271]. 

[4] M. Bruni, S. Matarrese, S. Mollerach and S. Sonego, Class. Quant. Grav. 
14 2585 (1997) [arXiv:gr-qc/9609040]; S. Matarrese, S. Mollerach and 
M. Bruni, Phys. Rev. D 58, 043504 (1998) [arXiv:astro-ph/9707278]; 
M. Bruni, F. C. Mena and R. K. Tavakol, Class. Quant. Grav. 19, L23 
(2002) [arXiv:gr-qc/0 107069]. 

[5] V. Acquaviva, N. Bartolo, S. Matarrese and A. Riotto, Nucl. Phys. B 
667, 119 (2003) [arXiv:astro-ph/0209156]. 

[6] J. M. Maldacena, JHEP 0305, 013 (2003) [arXiv:astro-ph/0210603]. 



42 



[7] N. Bartolo, E. Komatsu, S. Matarrese and A. Riotto, Phys. Rept. 402, 
103 (2004) [arXiv:astro-ph/0406398]. 

[8] G. Rigopoulos, Class. Quant. Grav. 21, 1737 (2004) [arXiv:astro- 
ph/0212141]. 

[9] K. A. Malik and D. Wands, Class. Quant. Grav. 21, L65 (2004) 
[arXiv:astro-ph/0307055]. 

[10] H. Noh and J. c. Hwang, Phys. Rev. D 69, 104011 (2004). 

[11] F. Vernizzi, Phys. Rev. D 71, 061301 (2005) [arXiv:astro-ph/0411463]. 

[12] K. Tomita, Phys. Rev. D 71, 083504 (2005) [arXiv:astro-ph/0501663]. 

[13] D. S. Salopek and J. R. Bond, Phys. Rev. D 42, 3936 (1990). 

[14] G. L. Comer, N. Deruelle, D. Langlois and J. Parry, Phys. Rev. D 49 
(1994) 2759. 

[15] N. Deruelle and D. Langlois, Phys. Rev. D 52, 2007 (1995) [arXiv:gr- 
qc/9411040]. 

[16] M. Sasaki and T. Tanaka, Prog. Theor. Phys. 99, 763 (1998) [arXiv:gr- 
qc/9801017]. 

[17] D. H. Lyth and D. Wands, Phys. Rev. D 68, 103515 (2003) [arXiv:astro- 
ph/0306498]. 

[18] G. I. Rigopoulos and E. P. S. Shellard, Phys. Rev. D 68, 123518 (2003) 
[arXiv:astro-ph/0306620] . 

[19] G. I. Rigopoulos and E. P. S. Shellard, JCAP 0510, 006 (2005) 
[arXiv:astro-ph/0405185]. 

[20] G. I. Rigopoulos, E. P. S. Shellard and B. W. van Tent, Phys. Rev. D 
73, 083521 (2006) [arXiv:astro-ph/0504508]. 

[21] G. I. Rigopoulos, E. P. S. Shellard and B. W. van Tent, Phys. Rev. D 
73, 083522 (2006) [arXiv:astro-ph/0506704]. 

[22] D. H. Lyth, K. A. Malik and M. Sasaki, JCAP 0505, 004 (2005) 
[arXiv:astro-ph/0411220]. 

[23] D. H. Lyth and Y. Rodriguez, Phys. Rev. D 71, 123508 (2005) 
[arXiv:astro-ph/0502578]; ibid., Phys. Rev. Lett. 95, 121302 (2005) 
[arXiv:astro-ph/0504045] . 

43 



[24] G. F. R. Ellis and M. Bruni, Phys. Rev. D 40, 1804 (1989). 
[25] S. W. Hawking, Astrophys. J. 145, 544 (1966). 

[26] M. Bruni, G. F. R. Ellis and P. K. S. Dunsby, Class. Quant. Grav. 9, 
921 (1992). 

[27] C. Gordon, D. Wands, B. A. Bassett and R. Maartens, Phys. Rev. D 
63, 023506 (2001) [arXiv:astro-ph/0009131]. 

[28] S. Groot Nibbelink and B. J. W. van Tent, arXiv:hep-ph/0011325; 
S. Groot Nibbelink and B. J. W. van Tent, Class. Quant. Grav. 19, 
613 (2002) [arXiv:hep-ph/0107272]. 

[29] D. Langlois, Phys. Rev. D 59, 123512 (1999) [arXiv:astro-ph/9906080]. 

[30] K. A. Malik, JCAP 0511, 005 (2005) [arXiv:astro-ph/0506532]. 

[31] R. M. Wald, "General Relativity", Chicago University Press, Usa 
(1984). 

[32] M. Sasaki and E. D. Stewart, Prog. Theor. Phys. 95, 71 (1996) 
[arXiv:astro-ph/9507001]. 

[33] A. A. Starobinsky, JETP Lett. 42, 152 (1985) [Pisma Zh. Eksp. Teor. 
Fiz. 42, 124 (1985)]. 

[34] M. Sasaki, Prog. Theor. Phys. 76, 1036 (1986); V. F. Mukhanov, 
Zh. Eksp. Teor. Fiz. 94, 1 (1988) [Sov. Phys. JETP 68, 1297 (1988)]. 

[35] A. Taruya and Y. Nambu, Phys. Lett. B428, 37 (1998). 

[36] J. M. Bardeen, Phys. Rev. D 22, 1882 (1980). 

[37] J. M. Bardeen, P. J. Steinhardt and M. S. Turner, Phys. Rev. D 28, 679 
(1983). 

[38] K. Enqvist, J. Hogdahl, S. Nurmi and F. Vernizzi, Phys. Rev. D 75, 
023515 (2007) [arXiv:gr-qc/0611020]. 

[39] F. Finelli, G. Marozzi, G. P. Vacca and G. Venturi, Phys. Rev. D 74, 
083522 (2006) [arXiv:gr-qc/0604081]. 

[40] K. Nakamura, Phys. Rev. D 74, 101301 (2006) [arXiv:gr-qc/0605107]; 
ibid., arXiv:gr-qc/0605108. 



44 



[41] K. Enqvist and A. Vaihkonen, JCAP 0409, 006 (2004) [arXiv:hep- 
ph/0405103]. 

[42] A. Jokinen and A. Mazumdar, JCAP 0604, 003 (2006) [arXiv:astro- 
ph/0512368]. 

[43] N. Barnaby and J. M. Cline, Phys. Rev. D 73, 106012 (2006) 
[arXiv:astro-ph/0601481]. 

[44] F. Vernizzi and D. Wands, JCAP 0605, 019 (2006) [arXiv:astro- 
ph/0603799]. 

[45] D. Seery, J. E. Lidsey and M. S. Sloth, arXiv:astro-ph/0610210. 



45 



